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ABSTRACT 

We present surface photometry and stellar kinematics of NGC 4431, a barred 
dwarf galaxy in the Virgo cluster undergoing a tidal interaction with one of its 
neighbors, NGC 4436. We measured its bar pattern speed using the Tremaine- 
Weinberg method, and derived the ratio of the corotation radius, D L , to the bar 
semi-major axis, as- We found D L /a B = 0.6±ol at 99% confidence level. Albeit 
with large uncertainty, the probability that the bar ends close to its corotation 
radius (i.e., 1.0 < D^jas < 1.4) is about twice as likely as that the bar is much 
shorter than corotation radius (i.e., Dl/cib > 1.4). 

Subject headings: galaxies: individual (NGC 4431) — galaxies: kinematics and 
dynamics — galaxies: elliptical and lenticular, cD — galaxies: photometry - 
galaxies: structure 

1. Introduction 

Stellar bars are observed in a substantial fraction of nearby disk galaxies (Knapen et al. 
2000; Eskridge et al. 2000; Marinova & Jogee 2006). Their growth is partly regulated by 
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the exchange of angular momentum with the stellar disk and the dark matter (DM) halo. 
For this reason the dynamical evolution of bars can be used to constrain the distribution 
of DM in the inner regions of galaxy disks. Using perturbation theory Weinberg (1985) 
predicted that a bar would lose angular momentum to a massive DM halo through dynamical 
friction, slowing down in the process. This prediction was subsequently confirmed in iV-body 
simulations (Debattista & Sellwood 1998, 2000; Athanassoula 2003; O'Neill & Dubinski 
2003; Sellwood & Debattista 2006). They found that bars are slowed efficiently within a few 
rotation periods if a substantial density of DM is present in the region of the bar. On the 
other hand, if the mass distribution is dominated by the stellar disk throughout the inner 
few disk scale-lengths, then the bar remains rapidly rotating (here defined as the bar ending 
just inside its corotation radius) for a long time. Thus the accurate measurement of the bar 
pattern speed, Q p , provides a way to discriminate whether the central regions of the host 
galaxy are dominated by baryons or by DM. 

Bar pattern speeds are most usefully parametrized by the distance-independent ratio 
7Z = Dl/clb, where Dl is the corotation radius and Ob the semi-major axis of the bar. Bars 
which end near corotation (1 < 1Z < 1.4) are termed fast, while shorter bars (1Z > 1.4) are 
said to be slow. The best way to determine the bar rotation parameter TZ is by measuring Q p 
with the model-independent method developed by Tremaine & Weinberg (1984, hereafter 
TW). To date this method has been used to measure f2 p of mainly massive (V c > 150 
km s" 1 ) early-type barred galaxies (see Corsini 2004 and references therein) and it has 
been successfully tested against iV-body simulations (Debattista 2003; O'Neill & Dubinski 
2003). All measured bars are compatible with being fast, implying that the inner regions of 
these galaxies are not dominated by DM. The only exception appears to be the slow bar in 
NGC 2915, a blue compact dwarf galaxy, for which the TW method was applied to H i radio 
synthesis data (Bureau et al. 1999). 

The minimum possible amount of DM in disk galaxies can be determined from the 
rotation curve, by scaling up the mass-to-light ratios of the luminous components to fit the 
central velocity gradient (see Bosma 1999 for a review). However, the lack of any obvious 
transition region from the bulge/disk-dominated to the halo-dominated part of rotation 
curves severely limits the ability to obtain a unique decomposition and derive the central DM 
distribution (e.g., Corsini et al. 1999). The structure of DM halos is more directly revealed 
in galaxies where the luminous components are thought to give a nearly zero contribution to 
the mass budget, because this makes mass modeling easier and the derived DM distribution 
less uncertain. For this reason dwarf galaxies were considered ideal targets for studying the 
properties of DM halos and for testing if they have a central cuspy power-law mass density 
distribution as predicted by cold dark matter (CDM) cosmology (Navarro et al. 1996, 1997; 
Moore et al. 1998, 1999). However, recent extensive studies have shown that in most of 
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them the contribution of the stellar disk can account for the inner rotation curve and DM 
halos with a constant-density core provide a better fit to than those with a density cusp 
(e.g., de Blok & Bosnia 2002; Swaters et al. 2003; Spekkens et al. 2005). 

Nonetheless, the interpretation of these results has been controversial because of possible 
observational systematic effects, such as the poor spatial resolution of radio maps (van den 
Bosch & Swaters 2001) and slit misplacement in optical spectroscopy (Swaters et al. 2003) 
as well as non-circular (Hayashi & Navarro 2006), and off-plane gas motions (Swaters et al. 
2003b; Valenzuela et al. 2005). 

Since measurements of pattern speeds with the TW method do not require high spatial 
resolution or slits passing exactly through the galaxy center and refers to the stellar com- 
ponent only, such a measurement in a dwarf disk galaxy can test unambiguously for the 
presence of a dense DM cusp. In the present Letter we attempt such a measurement for 
NGC 4431. 



2. NGC 4431 

NGC 4431 (VCC 1010) is a small (1'.7 x l'.l [de Vaucouleurs et al. 1991, hereafter RC3]) 
and faint (B^ = 13.74 [RC3]) nucleated dwarf galaxy (Binggeli et al. 1985) classified as 
dSBO/a by Barazza et al. (2002) following an unsharp mask analysis of the galaxy. The 
presence of a bar had earlier been missed by virtue of the galaxy's low luminosity and the 
relatively small angle between the bar's and disk's position angle (PA). Both the presence 
of a strong bar and trailing arms (which are clearly visible in Fig. 3 of Barazza et al. 
2002) and the rotationally-supported kinematics (Simien & Prugniel 2002; Pedraz et al. 
2002) are indicative of a genuine disk galaxy. The total absolute magnitude of the galaxy 
is Mg T = —17.41 corrected for inclination and extinction (RC3) and adopting a distance 
of 16.2 Mpc (Jerjen et al. 2004). It is located close to the center of the Virgo cluster at a 
projected distance of ~ 300 kpc from M 87. Its closest neighbor in projection is the dE6/dS0 
galaxy NGC 4436, at a projected separation of 3'.7 or ~ 18 kpc. 

3. Broad-band imaging 

Deep broad-band imaging of NGC 4431 was carried out with the Very Large Telescope 
(VLT) at the European Southern Observatory on April 1, 2000 as part of a project aimed at 
measuring distances, metallicities, and ages of dwarf galaxies in the Virgo cluster. Details 
of observations, data reduction, and analysis are reported elsewhere (Barazza et al. 2003; 
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Jerjen et al. 2004). 

We analyzed the R— band image of NGC 4431 obtained by Barazza et al. (2003), with 
foreground stars and background galaxies subtracted. We fitted ellipses to galaxy isophotes 
with the IRAF0 task ELLIPSE, after masking the northeastern part of the image to minimize 
the light contamination from NGC 4436, the outskirts of which fall in the field of view of 
the image. We first fitted ellipses allowing their centers to vary to test whether patchy 
dust obscuration and tidal deformations were present. Having found no evidence of varying 
ellipse centers within the errors of the fits for the inner 40" we concluded that in these regions 
there is little or uniform obscuration and no tidal deformations, and that the galaxy is a 
viable candidate for the TW method. At larger radii, the shift of ellipse centers (~ 1") of 
NGC 4431 in the direction of NGC 4436 may be due to the interaction between the two 
galaxies. The ellipse fits were then repeated with the ellipse center fixed and the resulting 
surface brightness profile is plotted in Fig. Id. The inclination (i = 48?7 ± 1?5) and disk 
PA (PAdisk = 6?9 ± 1?0) were determined by averaging the values measured between 50" and 
70". They are consistent within the errors with values obtained by fitting ellipses with free 
centers. 

We measured ag (Fig. 1) using four independent methods based on Fourier amplitudes 
(Aguerri et al. 2000), Fourier and ellipse phases (Debattista et al. 2002), and a parametric 
decomposition of the surface brightness profile (Prieto et al. 2001). The mean of the resulting 
values is our best estimate of a B and we assume the largest deviation from the mean as our 
error estimate (as = 21"9 ± 1"5). The bar length is in agreement with the radius where the 
spiral arms seen in the unsharp mask image of Barazza et al. (2002) start. 

4. Long-slit spectroscopy 

The spectroscopic observations of NGC 4431 were carried out in service mode at VLT 
on May 4-6, 2003 (run 1), April 21-22, 2004 (run 2) and May 18-23, 2004 (run 3). The Focal 
Reducer Low Dispersion Spectrograph 2 (FORS2) mounted the volume-phased holographic 
grism GRIS_1400V+18 with 1400 grooves mm -1 and the 0"7 x 6'.8 slit. The detector was 
a mosaic of 2 MIT/LL CCDs. Each CCD had 2048 x 2068 pixels of 15 x 15 /mi 2 . The 
wavelength range from 4560 to 5860 A was covered with a reciprocal dispersion of 0.645 A 
pixel -1 and a spatial scale of 0.250 arcsec pixel -1 after a 2 x 2 pixel binning. We obtained 
spectra with the slit at PA = 7?9 crossing the galaxy center (3 x 45 minutes in run 1) and 
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offset by 5"0 eastward (4 x 30 minutes in run 3). We also obtained spectra with the slit at 
PA = 6?9 and shifted by 5"0 eastward (2 x 40 minutes in run 2; 2 x 40 minutes in run 3) and 
westward (4x40 minutes in run 3) with respect to the galaxy nucleus. Two different PAs were 
chosen for the slits to test for the sensitivity of the TW measurement to the ~ 1° uncertainty 
in PAdisk- Comparison lamp exposures obtained for each observing night ensured accurate 
wavelength calibrations. Spectra of G and K giant stars served as kinematical templates. 
The average seeing FWHM was 1"2 in run 1, Cf!8 in run 2, and C'9 in run 3. Using standard 
MIDA£§] routines, all the spectra were bias subtracted, flat-field corrected, cleaned of cosmic 
rays, corrected for bad pixels, and wavelength calibrated as in Debattista et al. (2002). The 
accuracy of the wavelength rebinning 1 km s _1 ) was checked by measuring wavelengths 
of the brightest night-sky emission lines. The instrumental resolution was 1.45 ± 0.01 A 
(FWHM) corresponding to o inst = 35 km s _1 at 5170 A. The spectra obtained in the same 
run along the same position were co-added using the center of the stellar continuum as 
reference. In the resulting spectra the sky contribution was determined by interpolating 
along the outermost ~ 30" at the edges of the slit and then subtracted. 



5. Pattern speed of the bar 

To measure the pattern speed of the bar, O p , we used the TW method. For slits 
parallel to the disk major-axis, this relates the luminosity-weighted mean position, X, to the 
luminosity- weighted mean velocity V = XQ P sin i where i is the galaxy inclination. 

To measure V for each slit, we first collapsed the two-dimensional spectrum along its 
spatial direction between —90" and +90" in the wavelength range between 5220 and 5540 
A, obtaining a one- dimensional spectrum. The value of V was then derived by fitting the 
resulting spectrum with the convolution of the spectrum of the G0III star SA0119458 and a 
Gaussian line-of-sight velocity profile by means of the Fourier Correlation Quotient method 
(Bender 1990, hereafter FCQ) as done in Aguerri et al. (2003, hereafter ADC03). We 
estimated uncertainties by Monte Carlo simulations with photon, read-out and sky noise. 
As a precision check we analyzed separately the two eastward-offset spectra (Y = +5"0) 
which were obtained with the slit at PA = 6?9 in run 2 and 3, respectively. The two 
resulting estimates of V are in agreement within the errors (Fig. 2). To compute X for each 
slit, we extracted the luminosity profile from the R— band image along the position of the slit 
after convolving the image to the seeing of the spectrum. The R— band profile matches very 
well that obtained by collapsing the spectrum along the wavelength direction, supporting 
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that the slit was placed as intended. We used the R— band profiles to compute X because 
they are less noisy than those extracted from the spectra, particularly at large radii. By 
comparing the slits at PA = 6?9 and at PA = 7?9 we also verified that the error introduced 
by the ~ 1° uncertainty in disk PA is small for both X and V. We obtained Q p by fitting a 
straight line to the values of X and V to the data for slits at Y — ±5'/0 and PA = 6?9 and 
at Y = 0" and PA = 7?9 (Fig. 2). Although the Y = 0" slit was obtained at different PA, 
this slit constrains only the zero point (i.e., systemic) velocity and any slit passing through 
the center can be used. This gives f2 p sin 2 = 5.56 ± 1.39 km s _1 arcsec -1 (70.8 ± 17.7 km s -1 
kpc- 1 ). 

We used the FCQ to measure the line-of-sight velocity curve and velocity dispersion 
profile of the stellar component along the major axis (Fig. 3). Our measurements at PA = 6?9 
are in agreement within errors with those by Simien & Prugniel (2002) at PA = 6°. We 
derived the circular velocity in the disk region, V c — 94 ± 19 km s _1 , after a standard 
correction for the asymmetric drift as in ADC03. This value is in agreement with the circular 
velocity derived from the Tully-Fisher relation calculated in the R— band by Courteau (1997). 
Thus the corotation radius of the bar is Dl = V c /Vt v = l^!T^ 2 'l an d the ratio of the corotation 
radius to the bar semi-major axis K = D L /a B = 0.6tH- The error intervals on Dl and 1Z 
are 99% confidence level and were measured with Monte Carlo simulations as in ADC03. By 
excluding the possibility that DlI^b < 1.0 as this is thought to be unphysical (Contopoulos 
1980) we estimated that the probability that the bar ends close to its corotation radius (i.e., 
1.0 < D L /a B < 1.4) is 72%. This was calculated via Monte Carlo simulations by assuming a 
uniform distribution of a B and V c and Gaussian distribution of f2 p within their error ranges. 



6. Discussion 

By measuring the bar pattern speed of NGC 4431 we have demonstrated that the TW 
method is feasible using stars as tracer for dwarf galaxies. So far, this technique has been 
applied only to bright barred galaxies with the result that all the measured bars are consistent 
with being rapidly rotating (see Corsini 2004 for a review). For NGC 4431 we find a bar 
suggestive of being fast at 72% probability, although the measurement uncertainties in this 
galaxy preclude a stronger statement. 

As discussed by Debattista (2003) and Debattista & Williams (2004), the TW measure- 
ment of Q p is particularly sensitive to errors in the PA^. However, this is unlike to be 
an issue for NGC 4431. In fact, the values of X and V measured in the offset position at 
Y = +5'/0 from the spectrum at PA = 7?9 are in agreement within the errors with those 
from the two spectra at PA = 6?9. So, the fact that the value of the bar speed parameter 
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of NGC 4431 is nominally less than unity cannot be attributed to the uncertainty on PA disk 
but to the scatter in the values of X and V. 

A fast bar in NGC 4431 would suggest a common formation mechanism of the bar both 
in bright and dwarf galaxies. If disks were previously stabilized by massive DM halos, we 
exclude that these bars were produced by tidal interactions because they would be slowly 
rotating (Noguchi 1999). But, this is not the case even for galaxies, like NGC 4431 and the 
SBO NGC 1023 (Debattista et al. 2002) which show signs of weak interaction with a close 
companion, without being significantly perturbed. This finding allows us to conclude that 
there is no difference between TW measurements in isolated or mildly interacting barred 
galaxies (ADC03; Corsini et al. 2003). 

According to high-resolution iV-body simulations of bars in cosmologically motivated, 
CDM halos (see Debattista 2006 for a review) a fast bar would imply that NGC 4431 has 
a massive disk and does not reside inside a centrally-concentrated DM halo. This is in 
agreement with other studies of dwarf galaxies (Swaters 1999; de Blok et al. 2001; McGaugh 
et al. 2001; de Blok & Bosma 2002; Marchesini et al. 2002; Swaters et al. 2003; Spekkens 
et al. 2005) who analyzed the ionized-gas rotation curves of up to about 200 objects. By 
assuming minimal disks and spherical symmetry most of the data are better described by 
DM halos with an inner density core rather than the cusp predicted by CDM. 

However, there is no consensus whether these results based on gas dynamics represent an 
actual problem for the CDM paradigm or whether they could be reconciled with it by taking 
a variety of observational uncertainties (e.g., slit positioning, spatial resolution), analysis 
problems (e.g., binning and folding of the rotation curves), and modelling assumptions (e.g., 
non-circular and off-planar gas motions) into account (de Blok et al. 2003; Swaters et al. 
2003; Spekkens et al. 2005). Therefore, we argue that the ongoing debate will benefit from 
stellar dynamics to constrain the inner mass distribution of DM with the direct measurement 
of fl p in a sample of dwarf barred galaxies by means of integral-field spectroscopy. 
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Fig. 1. — Bar semi-major axis length of NGC 4431. (a) Phase angle of m = 2 Fourier 
component, (b) Bar/interbar intensity ratio, (c) PA of the deprojected isophotal ellipses. 
(d) Surface brightness decomposition in bulge (dashed line), disk (dotted line), and bar 
(dashed-dotted line). The vertical lines show the value of as obtained with each model. 
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Fig. 2. — Pattern speed of the bar in NGC 4431. The kinematic integrals V are plotted 
as a function of the photometric integrals X. The best fitting straight line has a slope 
fipsku = 5.56 ± 1.39 km s _1 arcsec -1 . Open circles correspond to slits at PA = 7?9 (Y = 0" 
and Y = +5'. / 0), while filled circles are for slits at PA = 6?9 (Y = ±5"0). 
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Fig. 3. — (a) The major-axis radial profile of the stellar line-of-sight velocity dispersion 
fitted with an exponential profile at r > 15" (solid line), (b) The major-axis radial profile 
of the stellar line-of-sight velocity (after subtracting the systemic velocity V sys = 948 ± 3 
km s _1 ) and the V c shu curve (solid line) with errors (dotted lines) obtained by applying 
the asymmetric drift for r > 15" as in ADC03. In (a) and (b) the measured profiles are 
folded around the center, with filled circles and asterisks referring to the N (receding) and S 
(approaching) sides, respectively. 



